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The conventional fluid description of multi-component plasma, supplemented by an appropriate 
equation of state for the macroparticle component, is used to evaluate the longitudinal sound velocity 
of Yukawa fluids. The obtained results are in very good agreement with those obtained earlier 
employing the quasi-localized charge approximation and molecular dynamics simulations in a rather 
broad parameter regime. Thus, a simple yet accurate tool to estimate the sound velocity across 
coupling regimes is proposed, which can be particularly helpful in estimating the dust-acoustic 
velocity in strongly coupled dusty (complex) plasmas. It is shown that, within the present approach, 
the sound velocity is completely determined by particle-particle correlations and the neutralizing 
medium (plasma), apart from providing screening of the Coulomb interaction, has no other effect 
on the sound propagation. The ratio of the actual sound velocity to its “ideal gas” (weak coupling) 
scale only weakly depends on the coupling strength in the fluid regime, but exhibits a pronounced 
decrease with the increase of the screening strength. The limitations of the present approach in 
applications to real complex plasmas are briefly discussed. 

PACS numbers: 52.27.Lw, 52.35.Dm 


I. INTRODUCTION 

Since the theoretical discussion of low-frequency dust- 
acoustic-waves (DAW) [H and their experimental discov¬ 
ery ii, the concept of dust-acoustic velocity became 
one of the central concepts in the physics of dusty (com¬ 
plex) plasmas. The observations of DAWs and measure¬ 
ments of dust-acoustic (sound) velocities have been used 
to obtain useful information about dusty plasmas systems 
under various conditions [Ml (we mainly concentrate 
on three-dimensional particle clouds here, which is also 
reflected in the reference list). In particular, experimen¬ 
tally measured DAW dispersion relations and sound ve¬ 
locities have been repeatedly used to estimate the electri¬ 
cal charge of particles in dusty plasmas under laboratory 
and microgravity conditions |16l420l| . employing various 
experimental techniques. For other aspects of DAWs and 
their significance for the held of dusty (complex) plasmas 
see, for instance. Refs. [2ll - l^ and references therein. 

Due to the large charge on the dust particles, com¬ 
plex plasmas often occur in a strongly coupled (liquid) 
state [23 - [^ . The effects of strong coupling were ne¬ 
glected in the original derivation of DAW [I|, although 
the dispersion relation of the dust waves can clearly be se¬ 
riously modified at strong coupling. Not surprisingly, the 
influence of strong coupling has been considered by many 
researchers, from the theoretical point of view, using dif¬ 
ferent approaches. These include quasi-localized charge 
approximation (QLCA) generalized hydrody¬ 

namics [3lj,[33, local field correction description [s^, and 
“multicomponent kinetic theory” . Molecular dynam¬ 
ics (MD) simulations have also been performed to obtain 


‘Also at Joint Institute for High Temperatures, Russian Academy 
of Sciences, Moscow, Russia 


wave dispersion relations in strongly coupled Yukawa flu¬ 
ids [351 l36l | . From experimental point of view, the topic 
has also received some attention [^ . l3^ . Overall, the 

effect of strong coupling on weave phenomena in com¬ 
plex plasmas remains a very important current research 
topic 

The purpose of this paper is to investigate to which 
extent the standard fluid description of multi-component 
plasma, supplemented by an appropriate equation of 
state for the particle component, can be used for quan¬ 
titative estimations of the longitudinal sound velocity 
in Yukawa fluids. It is known that the simplistic fluid 
treatment cannot reproduce the entire dispersion curve 
of longitudinal waves at strong coupling, and becomes 
particularly irrelevant in the regime of short wavelength. 
On the other hand, it can be expected to become more 
adequate for long wavelengths and low frequencies, i.e. 
in the regime of our present interest. There is also a 
direct evidence of the applicability of the conventional 
hydrodynamic description to sufficiently long waves in 
Yukawa s yst ems, which comes from numerical simula¬ 
tions |4^. In the simplest fluid formulation the ef¬ 
fects of strong coupling are entirely accounted for by an 
appropriate equation of state. Here we make use of the 
recently proposed simple and accurate practical expres¬ 
sions for the internal energy and pressure of Yukawa flu¬ 
ids, applicable across coupling regimes [4l|. This allows 
us to obtain quantitative results and to perform detailed 
comparison with the results from other approaches. 

In addition, we discuss another relevant issue, concern¬ 
ing the relation between the sound velocity in the conven¬ 
tional Yukawa system (charged particles interacting via 
Yukawa potential immersed into neutralizing medium) 
and that in a single component Yukawa system (an imag¬ 
inary one-component system of particles interacting via 
the model Yukawa potential without any neutralizing 
medium). This question, from somewhat different per- 
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spective, has been recently considered in Ref. 

The system considered in the following represents a 
collection of point-like highly charged particles in the 
neutralizing medium composed of plasma electrons and 
ions. The particles interact via the pairwise Yukawa (or 
Debye-Hiickel) interaction potential 


derivation of DAW dispersion relation in Ref. [I| . In this 
formulation electrons and ions provide equilibrium neu¬ 
tralizing medium and are described by 

— eui'S/cj) = TiVui, (2) 

eneVcti = TeS/rie, (3) 


V{r) = ((3^/r)exp(-r/AD), (1) 

where Q is the particle charge, Ad is the Debye screening 
length, and r is the distance between a pair of particles. 
Screening is produced by the neutralizing plasma and the 
screening length, in the linear approximation, is given by 
Ad = ADi/\/l + (TiUe/Teni), where Adi = y^Ti/dTre^rii 
is the ion Debye radius. Here ni(e) and T^(e) are the ion 
(electron) density and temperature, respectively, and e 
is the elementary charge. Since quasineutrality implies 
rii ^ rie (for negatively charged particles) and the elec¬ 
tron temperature is much higher than that of the ions 
(Te ^ Ti) in many practical situations, the screening is 
mostly associated with the ion component. 

It is conventional to characterize the state of Yukawa 
systems in terms of two dimensionless parameters. The 
first is the coupling parameter, T = Q^/oT, where 
a = (3/4^np)V3 is the Wigner-Seitz radius and T is 
the kinetic temperature of the particle component (tem¬ 
peratures are measured in energy units throughout the 
paper). This is roughly the ratio of the (bar Coulomb) 
interaction energy between neighboring particles to their 
kinetic energy. The second is the screening parameter 
K = a/ Ad, which is roughly the ratio of the interparticle 
separation to the screening length. 

Clearly, the idealized model described above over¬ 
simplifies considerably the actual rather complex in¬ 
teractions between the particles in real complex plas¬ 
mas [ 2 ^ 0,0 , although some experimentally observed 
trends can be reproduced by this simple consideration, 
at least qualitatively. However, more important in the 
present context is that simplifications involved make it 
possible to perform a direct comparison with the results 
obtained using other approaches (such as QLCA and 
MD). The limitations of our simplified model to describe 
real complex plasmas will be briefly discussed towards 
the end of the paper. 

The paper is organized as follows. In Section |lT] sim¬ 
plest fluid description of Yukawa model systems is for¬ 
mulated and an expression for the sound velocity is de¬ 
rived. Thermodynamic quantities of Yukawa fluids, nec¬ 
essary to evaluate the sound velocity, are provided in 
Section IHIl The obtained results for the sound velocity 
are then analyzed and benchmarked against previously 
published results in Section El This is followed by con¬ 
cluding remarks in Section fVl 


II. SIMPLEST FLUID DESCRIPTION 

We adopt the simplest fluid description of multi- 
component plasmas, similar to that used in the original 


where cj) is the electric potential. Equations (I2|) and © 
result in equilibrium Boltzmann relations for the ion and 
electron densities in the wave potential. 

The continuity and momentum equation for the parti¬ 
cle component are 

^+V(npvp) = 0, (4) 


avp , , VP 

^ + (vp • V)vp =-, 

ot rup rriprip 


(5) 


where rip and Vp is the particle density and velocity, nip 
is the particle mass, and P is the pressure associated with 
the particle component. 

In the limit of long-wavelength perturbations (acoustic 
regime) the system is quasineutral. 


en[ — eug -|- Qup = 0 . 


( 6 ) 


The standard linearization procedure is then applied 
to the set of equations ©-© . We assume nj = rijo + uji 
(j = e,i,p), (j) = (1)1, Vp = Vpi, and P = Pq + Pi, 
where the quantities with the subscript “1” correspond 
to small perturbations in the sound wave. All pertur¬ 
bations are proportional to exp(fkr — iujt), where lu is 
the wave frequency, and k is the wave vector. Since the 
sound wave is adiabatic [i^, the small change in pres¬ 
sure Pi is related to the small change in particle density 
npi by Pi = {dP/dnp)snpi, where S is the system en¬ 
tropy. The adiabatic compressibility can be expressed in 
terms of the isothermal compressibility using the ther¬ 
modynamic relation {dP/dnp)s = 'y{dP/dnp)T, where 
7 = C'pjCi/ is the adiabatic index. This results in the 
dispersion relation of the acoustic type 

, ,2 

2 _ 2 \ 2 I 2 

Cs = ■p' — ^p-^D + (7) 

where Cg is the sound velocity, ujp = y/AnQ^np/nip is 
the plasma-particle frequency, uxp = i/T/mp is the par¬ 
ticle thermal velocity, and ^ = {\/T){dP/dnp)T is the 
isothermal compressibility modulus. In the limit of cold 
particle system (T = 0) equation ([7]) reduces to the con¬ 
ventional dust-acoustic velocity 


Co = WpAD = 



1 Ht 

/ rUp \ 

1 1 + t + Ht' 


( 8 ) 


derived originally by Rao, Shukla and Yu [l| (see also 
[ 2 ^). Here Z = \Q/e\ is the particle charge number, 
r = TejTi is the electron-to-ion temperature ratio, and 











3 


H = Zn-plue is the so-called Havnes parameter. Note 
that normally the sound velocity Cg is much larger than 
the particle thermal velocity uxp due to a large factor 
Z. Thus, the assumption of cold particles is justified at 
weak coupling, since in this regime p, ^ 1 and 7 ~ 0 {\). 
(Note, however, that strictly speaking the assumption of 
weak coupling and the condition of cold particles are not 
consistent, since F —>■ 00 when T —>■ 0). The focus of 
the present study is mostly on the regime, when fj, can 
considerably deviate from unity. 

Rewritten in terms of reduced Yukawa state variables, 
K and r, the dispersion relation © becomes 

® 

This is identical to Eq. (27) from Ref. [i^ in the limit 
of vanishing wavenumbers (long wavelengths). The re¬ 
maining step to identify the influence of strong coupling 
on the sound velocity is to take the appropriate values 
for 7 and The proper choice is discussed in the next 
Section. 


III. THERMODYNAMIC FUNCTIONS OF 
YUKAWA FLUIDS 


Main thermodynamic quantities of interest here are 
the internal energy U, Helmholtz free energy E, pressure 
P, specific heats Cp and Cy, and the isothermal com¬ 
pressibility modulus fi = T~^{dP/dnp)x- If the internal 
energy is known, the following thermodynamic identities 
can be used to obtain other quantities |47l |: 


Cp 


U = -T^ 


d F 
ffTT 

\dV 


P = - 


-Cy = -T 


{dP/dT)l 

{dP/dV)T 


(lOa) 

( 10 b) 

( 10 c) 

(lOd) 


It is convenient to work with reduced units: u = U/NT^ 
f = F/NT, and p = PV/NT, Cp = Cp/N, Cv = Cy/N 
and to express all the derivatives in terms of Yukawa 
system phase state variables, k and T. To do that we 
assume that the unperturbed electron and ion densities 
are not related to the particle density and their tempera¬ 
tures are fixed. This implies T oc {aT)~^ oc T~^ny^ and 
K oc a oc rip This results in 

ar _ r sr _ i r _ n dn _ i n 

dT T ’ drip 3 rip ’ dT ’ dup 3 rip 

Thermodynamic identities yield the following re¬ 
lations between the reduced thermodynamic functions in 


(k, r) variables: 


f{K, T) = /id + r dr' [u(ir, r') - 3/2] /r', 
Jo 

rdf Kdf 


p{K,r) = 


3 dr 3dn 


( ll a) 

( ll b) 


where /id = In [(27r?i^/nipr)^/^np] — 1 is the ideal gas 
contribution to the free energy. The specific heat per 
particle at constant volume is expressed in terms of the 
reduced internal energy, 


cv(k, r) = u — r{du/dr). 


( 12 ) 


The isothermal compressibility modulus is related to the 
reduced pressure via 


/r(K,r) =p + 


r dp K dp 

3 (?r 3 dK 


(13) 


Finally, the difference between reduced specific heats at 
constant pressure and volume is 


Cp Cv — 


[p-r{dp/dryf 


(14) 


Note that it is conventional to decompose the thermo¬ 
dynamic quantities into the ideal gas contribution and 
excess contribution, associated with particle-particle cor¬ 
relations. We have not done this here in order to simplify 
the notation. The ideal gas contributions to the reduced 
energy and pressure are um = 3/2 and pid = 1, respec¬ 
tively. 

The multi-component character of the system under 
consideration implies that the thermodynamic quantities 
considered above can contain two distinct contributions. 
The first comes from particle-particle interactions and 
characterizes a single component Yukawa system without 
neutralizing background (it also contains the ideal gas 
contribution for non-interacting particles). The second is 
the plasma-related contribution, describing plasma and 
plasma-particle interactions. In particular, the plasma- 
related contribution to the internal energy of the system 
is M 


3r kF 


(15) 


The first term represents the (free) energy of the electron- 
ion plasma that, on average, neutralizes the charge of 
the particles, while the second term gives the (free) en¬ 
ergy of the plasma sheath around each particle, in the 
linear approximation. This latter term does not affect 
neither pressure, nor the compressibility modulus (and 
hence there is no effect on specific heats) [i^, It 
results, however, in the polarization force acting on the 
particles in non-uniform plasmas, see Appendix CT T he 
plasma-related contributions to p and p are thus j4l| 

Ppi = —3F/2it^ end /ipi = —3 F/k^. (16) 
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The necessity to account for the plasma-related contribu¬ 
tion can be easily understood realizing that it ensures the 
expected ideal gas values p = \ and /i = 1 in the high- 
temperature limit, when particle-particle correlations are 
completely absent. 

Thermodynamic properties of Yukawa systems have 
been extensively studied using various computational and 
analytical techniques. Relevant examples include Monte 
Carlo (MC) and molecular dynamics (MD) numerical 
simulations, as well as integral equation theoreti¬ 
cal studies [54l - l57j . Accurate values of various thermody¬ 
namic quantities have been tabulated in a wide region of 
(k, r) phase space (in particular, see Refs. [52|,[5^). Sim¬ 
ple practical expressions for u, p, and p have been put 
forward recently [dlj. These expressions stem from the 
original observation by Rosenfeld and Tarazona 
that the thermal component of the internal energy of var¬ 
ious soft repulsive systems (including the Yukawa case) 
exhibits quasi-universal dependence on the properly nor¬ 
malized coupling parameter. The resulting expressions 
are applicable in a wide range of coupling and demon¬ 
strate remarkable agreement with the results from nu¬ 
merical simulations. The expression for the e nerg y of the 
single component Yukawa fluid suggested in is 


,(K,r) — - -I- e -I- 


k{k -i- i)r 


(k -I- 1) -I- (k — l)e2 


■(5(r/r„,)2/5, 
(17) 

where r„i denotes the coupling parameter at the fluid- 
solid phase transition and the subscript “pp” means that 
only the contribution coming from particle-particle cor¬ 
relations (which also includes the ideal gas term) is con¬ 
sidered. The functional dependence rni(K) can be ap¬ 
proximated by [illMl 


rm(K) ^ 


172 exp(aK) 

1 + aK+ ’ 


(18) 


where the constant a = (47r/3)^/^ ~ 1.612 is the ra¬ 
tio of the mean interparticle distance A = to the 

Wigner-Seitz radius a. It has been also suggested [4l| to 
use 6 = 3.2 and e = —0.1 in Eq. (IT71) . 

The expressions for p, p, and 7 can be easily derived 
using Eqs. (HHl-dlHl), as has been done in Ref. [4l|. Note 
that when plasma-related contribution is properly ac¬ 
counted for, p and p tend to unity in the ideal gas limit 
(vanishing correlations) and become negative at suffi¬ 
ciently strong coupling. The adiabatic index 7 tends to 
5/3 in the ideal gas limit, exhibit a discontinuity at mod¬ 
erate coupling (when p = 0) and then rapidly approach 
the asymptote 7 = 1 at strong coupling. This is illus¬ 
trated in Eigure[T] where the adiabatic index 7 is plotted 
as a function of the coupling parameter T for four repre¬ 
sentative values of k. The discontinuity in 7 may deserve 
special attention, but has no impact on the problem un¬ 
der consideration, because the product jp remains finite. 
It decreases monotonously from the ideal gas value 5/3 to 
highly negative values when E increases. In this regime 
the sound velocity should decrease compared to its ideal 



FIG. 1: (Color online) Adiabatic index 7 = Cp/Cv as a func¬ 
tion of the coupling parameter F for four values of the screen¬ 
ing parameter: k = 1.0 (red solid curve), k = 2.0 (blue dashed 
curve), K = 3.0 (black dotted curve), and k = 4.0 (green solid 
curve). The discontinuity in 7 is the consequence of plasma- 
related contribution to the thermodynamic quantities. 


gas value (evaluated at the same k and T values), as be¬ 
comes immediately obvious from Eq. 


IV. SOUND VELOCITY 


Now we can evaluate the sound velocity of Yukawa 
fluids in a broad range of coupling. Eirst, we write 


7M = M + 


[p-V{dp/dV)f 
u-T{du/dT) ’ 


(19) 


where u, p, and p account for both particle-particle cor¬ 
relation and plasma-related effects, i.e. u = Upp -I- Upi, 
p = Ppp -I- Ppi, and p = Ppp + Ppi- The first important 
observation is that when substituting Eq. HUD into the 
expression for the sound velocity ®, the plasma-related 
contribution to the isothermal compressibility modulus 
Ppi = —3r/K^ cancels out exactly the plasma contribu¬ 
tion to the dispersion relation given by the term 1 /k^ in 
®. The expression for the sound velocity can be rewrit¬ 
ten as 



[p-ndp/dT)f 

3r[u-r(5M/5r)]y 


( 20 ) 


The next important observation is that since the plasma- 
related contributions Upi and Ppi are both linear in T, they 
have no effect on the second term in brackets of Eq. O- 
Thus, the sound velocity of a system of charged particles 
immersed in the neutralizing plasma environment is equal 
to that of an imaginary single component Yukawa system. 
In other words, the magnitude of the sound velocity is 
completely determined by particle-particle correlations 
and the neutralizing medium only affects the interpar¬ 
ticle interactions, but has no other effect on the sound 
propagation. Therefore, the expressions for Upp(K,r), 
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FIG. 2: (Color online) Reduced sound velocity of Yukawa 
fluids, Cs/cUpfl, as a function of the screening parameter k. 
The solid curves correspond to the results of the simple fluid 
approach of this paper for F = 10 (blue curve) and F = 100 
(red curve). The dashed curve is plotted using QLCA result 
of Ref. given by Eqs. (I21II and (1221) . 


Ppp(K,r), and /rpp(K,r) proposed in Ref. [4l| can be di¬ 
rectly substituted into Eq. dini) to evaluate the sound 
velocity for a given pair of k and F. We have done this, 
the obtained results are discussed below. 

First, it is convenient to benchmark our results against 
quantitative data published in previous works. For exam¬ 
ple, Kalman et al. suggested the following expression 
for the longitudinal sound velocity, based on the results 
of QLCA model, 

Cs = WpO [1 /k^- b/( k)] , (21) 


where 


/(k) ~ -0.0799 - 0.0046^2 -b O.OOIGk'^ (22) 

is a fitting function applicable for k < 2.5 (note, that 
in this approach Cg/wpO is a function of the screening 
parameter k alone). In Figure [2] we plot the correspond¬ 
ing curve along with our calculation for the two values 
of the coupling parameter, F = 10 and F = 100. The 
curve corresponding to stronger coupling is practically 
indistinguishable from the QLCA result of Ref. [HU . 

A more detailed comparison with the QLCA is pro¬ 
vided in Table m The first two columns specify the loca¬ 
tion of the system in terms of screening parameter k and 
reduced coupling parameter F/Fm. The third column 
lists the values of the reduced sound velocity obtained 
using the QLCA model in Ref. [2^. The last column 
contains the values of the reduced sound velocity calcu¬ 
lated from the simple fluid approach of this paper us¬ 
ing Eq. (l20l) . Some weak dependence on the coupling 
strength is now present. The sound velocity obtained 
using the simple fluid approach is slightly smaller than 
that from the QLCA method, but the overall agreement 
is quite good. It is unlikely that the difference arises 
due to an approximate character of the equation of state 
employed here. The latter agrees with the accurate nu¬ 
merical data to within a tiny fraction of a percent in the 


TABLE I: Reduced sound velocity Cs/utpa of Yukawa fluids 
as calculated from the QLC approximation and present fluid 
model for several phase state points. QLCA data are from 
Ref. [^. For details see the text. 


K 

r/r,n 

QLCA 

Fluid 

1.0 

0.12 

0.96 

0.95 

1.0 

0.70 

0.96 

0.94 

2.0 

0.12 

0.42 

0.41 

2.0 

0.70 

0.41 

0.39 

3.0 

0.12 

0.23 

0.21 

3.0 

0.70 

0.21 

0.19 



q 


FIG. 3: (Color online) Long-wavelength dispersion of the lon¬ 
gitudinal waves in Yukawa fluids near freezing, plotted in the 
{q,ijj/u!p) plane, where q — ka is the reduced wavenumber. 
Symbols correspond to the results from numerical simulations 
of Ref. [ 3 ^ for k = 1.0 (red circles), k = 2.0 (blue triangles), 
and K = 3.0 (olive rhombuses). The corresponding dashed 
lines correspond to the acoustic asymptotes uj = kcs with the 
sound velocity Ca calculated using the fluid approach described 
in the present paper [Eq. (1201) ]. 


regime k 5 . 3 and F/Fm > 0.1 [4l|. Note that the QLCA 
model |62| is by construction more appropriate to de¬ 
scribe high-frequency short-wavelength phenomena and 
as such it is not necessarily more accurate than the fluid 
approach in predicting the sound velocity. 

In Figure[3]we compare the long-wavelength part of the 
longitudinal dispersion relation of Yukawa fluids obtained 
in a numerical (MD) experiment in Ref. (35j | with the 
acoustic asymptote uj = kcs, where Cs is calculated from 
the present fluid approach. The behavior of the numeri¬ 
cally obtained dispersion curves is clearly consistent with 
the respective acoustic asymptotes (dashed lines). More¬ 
over, Figure [3] demonstrates that at sufficiently strong 
screening (say k > 2 ), the acoustic asymptote describes 
well the dispersion curve up to g ~ 1, where q = ka is 
the reduced wavenumber. This would be sufficient for 
many experimental investigations of DAWs in complex 
plasmas. 

Having benchmarked the present fluid results against 
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K 

FIG. 4: (Color online) Contour plot of the reduced sound 
velocity Cs/co of Yukawa fluids in the plane (fi:,r/rin). Cal¬ 
culations are made using the fluid model described in this 
paper. 

previous results from the QLCA model and MD simu¬ 
lations, let us investigate the dependence of the sound 
velocity on coupling and screening in detail. It is par¬ 
ticularly useful to analyze the behavior of the quantity 
Cs/co, which is the ratio of the actual sound velocity of a 
Yukawa fluid to the respective limiting “ideal gas” (weak 
coupling limit) value given by Eq. (| 8 ]). The contour plot 
of this quantity in the (k, E/r„i) plane is shown in Fig. 2] 
First, we observe only weak dependence of the quantity 
Cs/cq on r deep in the fluid regime (on approaching the 
fluid-solid phase transition). This implies that the ab¬ 
solute value of Cs increases with T, because cq oc v^. 
Second, we observe that the ratio Cs/cq is sensitive to 
the screening parameter, and decreases as k increases. 
It drops by almost one order of magnitude on the way 
from the weakly screened regime k ^ 1 to the strongly 
screened regime with k ~ 5. 

Figure [5] shows the dependence of Cs/cq on the reduced 
coupling parameter for four values of the screening pa¬ 
rameter (from K = 1.0 to K = 4.0). It demonstrates that 
the main drop in the reduced sound velocity occurs al¬ 
ready in the regime of relatively weak coupling. In the 
strong coupling regime (say F/Fm > 0.1) the decrease 
is very slow, which in some sense justifies neglecting the 
F-dependence, like in Eq. (EU above. 

V. CONCLUDING REMARKS 

The main results obtained in this study can be summa¬ 
rized as follows. We have applied the standard fluid de¬ 
scription of multi-component plasmas, supplemented by 
an appropriate equation of state, to evaluate the sound 



r/r„ 


FIG. 5: (Color online) The reduced sound velocity Cs/co of 
Yukawa fluids versus the reduced coupling parameter F/Fm. 
The curves (from top to bottom) correspond to re = 1.0, re = 
2.0, re = 3.0, and re = 4.0, respectively. Calculations are made 
using the fluid model described in this paper. The equation 
of state used in these calculations may not be very accurate 
at F/Fm < 0.1. 


velocity in Yukawa fluids. It turns out that the sound ve¬ 
locity is completely determined by particle-particle cor¬ 
relations. The obtained values of sound velocities are 
in rather good agreement with the previously published 
results obtained using QLCA approach and MD sim¬ 
ulations. The main trends observed include slow de¬ 
crease of the sound velocity of a Yukawa fluid compared 
to its “ideal gas” scale with increase in the coupling 
strength, and more pronounced decrease with increase 
of the screening strength. Overall, the standard fluid de¬ 
scription, with a proper model for an equation of state, 
provides (perhaps not very surprisingly) simple yet accu¬ 
rate tool to evaluate sound velocity in Yukawa fluids and 
related systems in a rather broad parameter regime. 

Among the possible applications of the present re¬ 
sults, the topics related to low-frequency wave phenom¬ 
ena in complex (dusty) plasmas seem particularly rele¬ 
vant. However, it is important to take into account the 
following circumstances, which can limit the applicability 
of the simple model discussed here in some practical situ¬ 
ations. First of all, when deriving the dispersion relation 
in Section ini the most simple possible formulation of the 
problem has been employed. In real situations one may 
need to account for the presence of external electric fields 
and resulting drifts between different charged species, 
various kinds of collisions present in the system, effects of 
plasma production and loss, particle charge variations in 
the presence of the waves, additional forces acting on the 
particle component, etc. Most of these effects can be rel¬ 
atively easily included into the conventional fluid formal¬ 
ism, for some relevant examples we refer to Refs. [63l - [70l| . 
Thus, careful analysis of the most important processes in 
each concrete practical situation is required. Inclusion of 
these processes into consideration should not become a 
major problem from the theoretical point of view. 

The second class of problems is related to the openness 
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of the complex plasma systems. Plasma electrons and 
ions are continuously lost on the particle surface and the 
particle charge is set by the condition of no net electrical 
current to the surface (or, equivalently, floating poten¬ 
tial at the particle surface). This is known to result in 
some deviations from the Yukawa-type potential around 
the particles [12, and, therefore, some deviations from 
the thermodynamic functions of conventional Yukawa flu¬ 
ids can also be expected. Perhaps even more important 
is that the particle charge in complex plasmas is not 
fixed, but depends on the parameters of the surround¬ 
ing plasma. In particular, the charge becomes a function 
of the particle density via the so called “charge canni¬ 
balism” effect This effect operates as follows: 

When the particle density increases, the negative charge 
carried by the particle component also increases, which 
results in some reduction of the electron-to-ion density 
ratio (electron depletion) in view of the quasineutrality 
condition. In turn, this suppress the efficiency of electron 
collection by the particle surface compared to that of the 
ions. The particle charge becomes less negative, i.e. de¬ 
creases in the absolute magnitude compared to the case of 
an individual particle. In general, the relation between 
the particle charge and number density and the densi¬ 
ties of electrons and ions in complex plasmas is governed 
by the quasineutrality condition and the competition be¬ 
tween specific plasma production and loss mechanisms 
operating in a given situation. All this indicates that the 
consideration of an idealized Yukawa system with fixed 
particle charges and background plasma density can be 
in many cases insufficient to mimic the actual thermody¬ 
namics of real complex plasmas. How large modihcations 
can be and whether they can be evaluated using conven¬ 
tional thermodynamic approaches require special careful 
investigation. We leave this for future work. 


are listed in Table HI We also thank Satoshi Hamaguchi 
for providing numerical data on wave dispersion relations 
in Yukawas fluids, partly reproduced in Figure [3] This 
study was partially supported by the Russian Science 
Foundation, Project No. 14-12-01235. 


Appendix A: Polarization force on an individual 
particle 


In the case of an individual particle, the only contribu¬ 
tion to the energy is that from the sheath around the par¬ 
ticle. The particle energy is therefore U = —T{kT/2) = 
—Q^/2Xd. If the plasma is non-uniform such a particle 
will be acted upon by the force 


Fpoi = -V[/ = -^^, (Al) 

which is known as the polarization force [Zllll- The 
polarization force is small in most practical situations oc¬ 
curring in complex (dusty) plasmas, but can affect con¬ 
siderably the dispersion of dust acoustic waves as has 
been pointed out in Ref. [62 |. 

In the present consideration we treated the particle 
charge Q constant. If this limitation is relaxed, then 
another contribution to the force appears formally 


^ QVQ 
Ad 


(A2) 


Acknowledgments 

We thank Zoltan Donko for providing the values of 
sound velocities obtained using the QLCA model, which 


which would push a particle in the region where its charge 
is higher. The question whether this is a real force or an 
artifact of this consideration and, in the first case, its sig¬ 
nificance for complex plasmas deserves separate detailed 
consideration. 
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